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Scalar-tensor theories of gravity provide a consistent framework to accommodate an ultra-light 
quintessence scalar field. While the equivalence principle is respected by construction, deviations 
from General Relativity and standard cosmology may show up at nucleosynthesis, CMB, and solar 
system tests of gravity. After imposing all the bounds coming from these observations, we consider 
the expansion rate of the universe at WIMP decoupling, showing that it can lead to an enhancement 
of the dark matter relic density up to few orders of magnitude with respect to the standard case. 
This effect can have an impact on supersymmetric candidates for dark matter. 
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I. INTRODUCTION 

According to our current understanding [l| , Dark Mat- 
ter (DM) and Dark Energy (DE) represent the two ma- 
jor components of the present Universe. Surprisingly, it 
is found that the DM and DE energy densities, poM and 
Pde, are today roughly the same (differing only by a fac- 
tor of two) , while their ratio has been varying by several 
orders of magnitude in the past history of the Universe. 

It seems quite natural, then, to explore the possibil- 
ity of a DM-DE interaction which could account for this 
coincidence. This approach, however, is not free from 
problems if the DE component is interpreted in terms of 
a dynamical quintessence scalar field. Indeed, such a 
scalar is constrained to be extremely light in order to fit 
the data, giving rise to unwanted long-range forces which 
may represent a severe threat to the equivalence princi- 
ple. In addition, couplings of the quintessence scalar with 
the gauge field strengths are potential sources of danger- 
ous time variations of the fundamental constants (30l|. 

A possible way of coping with a very light scalar while 
avoiding these shortcomings, is to choose to work in 
the framework of a scalar-tensor gravity (ST) theory [||, 
where by construction matter has a purely metric cou- 
pling with gravity. It has been shown 9] that in this case, 
the scalar field can benefit from an attraction mechanism 
which, during the matter dominated era, makes ST over- 
lapping with standard General Relativity (GR). At the 
same time, ST may possess a second "attraction mecha- 
nism" "9^ which will ensure the correct evolution oi pde 
along a so-called 'tracking' solution. 

While ST can very closely reproduce standard GR at 
the present time, it may however lead to major differ- 
ences in the past evolution of the Universe, differences 
which may result in observable consequences for us to- 
day. For example, it has been shown [ill IT^ 
that ST theories may have a profund impact on nucle- 
osynthesis. At the same time, a curious fact has recently 
come to attention: a non-conventional dynamics of the 
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quintessence scalar in the past history of the universe 
may remain 'hidden' to the available cosmological obser- 
vations, but manifest itself through the DM relic abun- 
dance [13, [la [13 • It is then worth studying if ST may 
provide a viable quintessence candidate and at the same 
time have an impact on the DM relic abundance. We 
have considered this possibility and computed explicitly 
the differences from the standard case. 

When considering ST theories, the departures from 
standard cosmology are mainly due to the different ex- 
pansion rate {H) which they determine. Such deviation 
from the usual expansion rate of GR bears potentially 
relevant consequences in all those phenomena which are 
closely dependent on the timing in which they occur. The 
aim of the present work is to study the possible modifica- 
tions of the expansion rate of the Universe in ST at the 
time of Cold Dark Matter (CDM) freeze-out, focussing 
on Weekly Interacting Massive Particles (WIMPs) as the 
most natural candidates for CDM. As it is well known, 
their present relic density depends on the precise moment 
they decouple and, in turn, on the precise moment the 
WIMP annihilation rate equals the expansion rate of the 
Universe. We expect then that a variation of H in the 
past may lead to measurable consequences on the WIMP 
relic abundance. 

In order to assess the allowed departure of H from its 
standard value at WIMPs freeze-out, we have to take 
into account the bounds imposed on ST by phenomena 
at later epochs We already mentioned the crucial 

test of nucleosynthesis: passing the BBN exam will im- 
ply rather severe restrictions on the coupling of the scalar 
field of ST to matter and also on its initial conditions at 
temperatures much higher than the WIMP freeze-out. 
Coming to later epochs, we have to consider the photon 
decoupling and the consequent restrictions imposed by 
CMB data Interestingly enough, we will show that 
the BBN "filter" on ST is so efficient that it drastically 
limits any visible effect on the CMB spectrum (in partic- 
ular, shifts in the peaks positions are forced to be smaller 
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than the experimental error). More restrictive than the 
CMB probe turn out to be the GR tests, in particular 
after the tight bound on the post-newtonian parameter 
7, recently provided by the Cassini spacecraft ^2^. This 
constraint becomes quite relevant when combined with 
that on the scalar equation of state of coming from 
SNe la data We will find that significant devia- 

tions from standard cosmology are possible only if the 
DE equation of state differs appreciably from —1 today. 
In other words, if DE is a cosmological constant, it is un- 
likely that future cosmological observations will be able 
to discriminate between ST and GR. 

The question we intend to explicitly tackle is the fol- 
lowing: taking all the abovementioned restrictions (BBN, 
CMB, GR tests) into account, how much can the Hubble 
parameter H at the time of WIMP freeze-out differ from 
its canonical value if ST replaces GR? In other words, 
how much is the WIMP relic density allowed to vary, if 
we consider ST instead of GR? 

We find that in ST theories the expansion rate of the 
Universe at few GcVs can profoudly differ from the usual 
value obtained in GR (with variations up to five orders of 
magnitude) and, yet, allow the correct light elements pro- 
duction at BBN. This situation is perfectly analogous to 
the 'kination' effect studied in 15], where a modification 
of H at WIMP freeze-out was induced by a short pe- 
riod of dominance of the scalar kinetic energy, although 
with some deal of fine-tuning. In the case considered 
here, the effect of ST on H depends on the strength of 
the scalar-matter coupling, however no particular fine- 
tuning is needed to pass the severe nucleosynthesis test 
even when large modifications of H at freeze-out occur. 
This means that the attraction of ST towards GR pro- 
ceeds very rapidly during the cooling of the Universe from 
the few GeVs of WIMPs freeze-out down to the MeV 
range of nucleosynthesis. The overlap of ST with GR 
can subsequently be very efficient leading to ST scenar- 
ios which can hardly be disentangled from ordinary GR 
in present tests at the post-newtonian level. The fact 
that ST strongly affects the number of CDM particles 
may turn out to be the major signature of these theories. 

^From the point of view of particle physics model 
building, these large variations in the WIMPs number 
density today is of utmost relevance. Particles which 
were not considered suitable to play a significant role 
in CDM scenarios can be rescued because of their en- 
hanced number density. On the other hand, particles (or 
regions of the parameter space for certain WIMPs can- 
didates), which in usual GR scenarios constitute promis- 
ing CDM candidates, would be excluded because their 
boosted number would overdose the Universe. These 
considerations become of particular interest if we focus 
on the case where the WIMPs correspond to the lightest 
supersymmetric particle. A complete analysis of the cos- 
mologically excluded and cosmologically interesting re- 
gions of the SUSY parameter spaces in different SUSY 
contexts, when ST is considered, is presently in progress 



II. SCALAR-TENSOR THEORIES OF GRAVITY 

ST theories represent a natural framework in which 
massless scalars may appear in the gravitational sector 
of the theory without being phenomenologically danger- 
ous. In these theories a metric coupling of matter with 
the scalar field is assumed, thus ensuring the equivalence 
principle and the constancy of all non-gravitational cou- 
pling constants |2^. Moreover, as discussed in a 
large class of these models exhibit an attractor mecha- 
nism towards GR, that is, the expansion of the Universe 
during the matter dominated era tends to drive the scalar 
fields toward a state where the theory becomes indistin- 
guishable from GR. 

ST theories of gravity are defined by the action 
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The matter fields are coupled only to the metric ten- 
sor and not to i.e. Sm = Sm,[^m, 9fj.u]- R is the 
Ricci scalar constructed from the physical metric g^u- 
Each ST model is identified by the two functions w(<I>) 
and U($). For instance, the well-known Jordan-Fierz- 
Brans-Dicke (JFBD) theory Q corresponds to u}{(^) — u 
(constant) and U(<i>) = 0. 

The matter energy-momentum tensor is conserved, 
masses and non-gravitational couplings are time inde- 
pendent, and in a locally inertial frame non gravitational 
physics laws take their usual form. Thus, the 'Jordan' 
frame variables g^v and <& are also denoted as the 'phys- 
ical' ones in the literature. On the other hand, the equa- 
tions of motion are rather cumbersome in this frame, as 
they mix spin-2 and spin-0 excitations. A more conve- 
nient formulation of the theory is obtained by defining 
two new gravitational field variables, and the dimen- 
sionless field (p, by means of the conformal transformation 
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Imposing the condition 
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the gravitational action in the 'Einstein frame' reads 
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and matter couples to ip only through a purely metric 
coupling, 



(6) 



In this frame masses and non-gravitational coupling con- 
stants are field-dependent, and the energy-momentum 
tensor of matter fields is not conserved separately, but 
only when summed with the scalar field one. On the 
other hand, the Einstein frame Planck mass is time- 
independent and the field equations have the simple form 
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where 



rpif 



and Tf^i, — 2{~g)^^/^ dSm/Sg'^'^ is the matter energy- 
momentum tensor in the Einstein frame. The rele- 
vant point about the scalar field equation in JT)) is that 
its source is given by the trace of the matter energy- 
momentum tensor, T = g^'^T^^, which implies the (weak) 
equivalence principle. Moreover, when a{ip) = the 
scalar field is decoupled from ordinary matter and the 
ST theory is indistinguishable from ordinary GR. 

We next consider an homogeneous cosmological space- 
time 



where the matter energy-momentum tensor admits the 
perfect-fluid representation 

T^'' = {p + p) u^'u" -vg^'" , 

with g^y u^u^ = 1. 

The Friedmann-Robertson- Walker (FRW) equations 
then take the form 
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with the Bianchi identity 

d{pa^) +pda^ ^ {p-3 p) a^d\ogA{p). 
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The physical proper time, scale factor, energy, and 
pressure, are related to their Einstein frame counterparts 
by the relations 



df = A{ip)dT, d = A{(p)a, p = A{(py^p, p = A{(py^ 
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Defining new dimensionless variables 
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and setting k = (flat space) the field equation of motion 
takes the more convenient form 
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where primes denote derivation with respect to N. This 
will be our master equation. 

The effect of the early presence of a scalar field on the 
physical processes will come through the Jordan-frame 
Hubble parameter H = dloga/dr: 



H = H 
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where H = dloga/dr is the Einstein frame Hubble pa- 
rameter. In the flat-space case (A; = 0), Eq H13|l finally 
gives: 
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III. EVOLUTION OF THE FIELD 

A. Radiation domination 

During radiation domination the scalar field energy 
density is suppressed, A <C 1, and the first term in the 
RHS of Eq. l(T^ is proportional to 
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where the sum runs over all particles in thermal equilib- 
rium, while pm is the contribution from the decoupled 
and pressureless matter abundance. 

During radiation domination, ptot — 7r^T*/30, where 
T is the Jordan-frame temperature, and the contribution 
from a single particle in equilibrium gives 
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with yA = niA/T, gA the number of degrees of freedom 
of A, g^^ the number of relativistic degrees of freedom and 
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FIG. 1: The function y'^F[y], with y = m/T defined by 
Eq. 1171 . The upper (lower) curve corresponds to bosons 
(fermions) . 



where ea = {va + x^Y^"^ and the minus (plus) sign 
in the denominator of the integrand holds for bosons 
(fermions). In Fig. we plot y^F[j/]. We see that it 
is different from zero only around y — 1, that is, for 
T ~ mA- For higher temperatures it is quadratically sup- 
pressed in y, approaching the relativistic regime in which 
1 — Swyi = 0. For lower temperatures it is Boltzmann- 
suppressed. Then, as emphasized in 0,0], the field Lp 
evolves even during radiation domination, receiving a 
'kick' each time a particle in equilibrium becomes non- 
relativistic. The second term in Eq. IjlSI) is suppressed 
as Toq/T, so it becomes relevant only as equivalence is 
approached. 

In the following, we will consider the evolution of the 
field from the freeze out of the DM particles down to 
today, so we will have to take into account all particles 
of masses between the freeze-out and matter-radiation 
equivalence. 



B. Matter domination 

During matter domination 1 — iw ~ 1 and, as long as 
the field energy density is subdoniinant (A ^ 1), the RHS 
of the equation of motion H12|l is given by —V2a{ip) and 
the field evolution depends on the form of the coupling 
function a{ip). As already mentioned, the JBD theory is 
given by a constant a, and the value a = corresponds 
to GR. A very attractive class of models is that in which 
the fimction a{(f) has a zero with a positive slope, since 
this point, corresponding to GR, is an attractive fixed 
point for the field equation of motion Q . 

It was emphasized in Ref. !| (see also that the 
fixed point starts to be effective around matter-radiation 
equivalence, and that it governs the field evolution until 
recent epochs, when the quintessence potential becomes 
dominant. If the latter has a run-away behavior, the 
same should be true for a{ip), so that the late-time be- 
havior converges to GR. 



C. Late-time beahavior 

The evolution of the field during the last redshifts de- 
pends on the nature of DE. We will consider two possi- 
bilities: a cosmological constant and a inverse-power law 
scalar potential for Lp, which can be collectively repre- 
sented by the potential 
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S — corresponding to the cosmological constant. 

In general, a cosmological constant in the Einstein 
frame does not correspond to a cosmological constant 
in the Jordan frame, as one can read from Eq. Q . How- 
ever, present tests of GR (see next section) imply that at 
late times A{(p) ~ 1, so that the two frames are almost 
coincident and the expansion histories during the last few 
redshifts are practically indistinguishable. 

For the purpose of this paper, that is the analysis of the 
impact of DE on ST cosmology, the situation in which the 
quintessence field is different from ip and decoupled from 
it would be basically the same as that with a cosmological 
constant, since in both cases the second term in the RHS 
of Eq. ifT^ vanishes. 

The main feature of the potential in (|18|l for (5 > is 
the existence of 'tracker' solutions, which are attractors 
in field space I'^ the a ^ limit, the late-time be- 
havior is completely determined by the two parameters 
A = A^/p^ and 5. A non- vanishing a would modify the 
behavior of the field today, hopefully keeping the desir- 
able property of insensitivity to the initial conditions. 

In this paper, we will consider the following choice for 



A{ip) = 1 + Be 



corresponding to 
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which has a run-away behavior with positive slope, as re- 
quired by the discussion at the end of the previous sub- 
section. The choice for the parameters B and /3 will be 
discussed in the following section. 

In Fig.|21we show the evolution of the background and 
field energy density. We see that field energy densities 
corresponding to different initial conditions converge to 
the same solution, driven by a{ip). Notice also that the a- 
attractor becomes effective even before matter-radiation 
equivalence, due to the non-rclativistic decoupling effect 
explained above. 

In Fig. ^ we show the region of parameter plane X-6 
giving w^p < —0.7, where 
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and 0.65 < < 0.75. We see that in the ST case, 
{B 7^ 0) the region giving more negative values of the 
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equation of state is somehow enlarged with respect to 
pure GR quintessence. However, in the observationally 
allowed region for fi^ the influence of the B parameter 
is negligible. 
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one, in which the Planck mass at nucleosynthesis was the 
same as today. We obtain: 



- ^GR 



-"gr 
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The change in the expansion rate is completely anal- 
ogous to that obtained by adding extra neutrinos to the 
standard GR case. Using the bound AA^ < 1 (which is 
more restrictive than those obtained for instance in i25j). 
we get 



A{ipo) 



< 1.08. 



(24) 



B. General relativity tests 



FIG. 2: Evolution of the energy density of the background 
(upper solid line) and of three typical solutions for the scalar 
field. We see that different initial conditions converge to the 
same solution. 
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FIG. 3: The regions in the X-S parameter plane giving Wip < 
-0.7 (dark grey) and 0.65 < < 0.75 (light grey). The left 
plot is the pure GR case (5 = 0) while the right one is for 
ST with B = 0.1, /? = 8. 



IV. PHENOMENOLOGICAL BOUNDS 

A. Nucleosynthesis 

Assuming aip' ~ in Eq. (|14|l ~ as we have checked 
numerically for the solutions relevant in this analysis ~ 
the Jordan frame expansion rate during nucleosynthesis 
may be approximated as 

ij2 ~ A2(09)^p. (22) 

The above expression should be compared to the GR 



At the post-newtonian level, the deviations from GR 
may be parametrized in terms of an effective field- 
dependent newtonian constant jsij 

G^G{^) = G,A{ipf{\ + a\ip)), 

and two dimensionless parameters ^pn and /?pjv which, 
in the present theories turn out to be [2^ 

2 9 

- 1 = , /3p^-l = ^__^, (25) 

where k — da/d^p. 

A new constraint on the parameter ^pn has been ob- 
tained recently using radio links with the Cassini space- 
craft [13, 

7P^-1:= (2.1±2.3) X 10"^ (26) 

Present bounds on (3pM — 1 are 0(10^"*) and are less 
restrictive for our choice of a(<^), since Koa^((po) — 

The bound from the Cassini spacecraft turns out to 
be quite strong when used in connection with the one on 
the equation of state w^p from SNe la. In Fig. Q we 
show the excluded region in the fi-Wtp plane implied by 
Eq. We see that an equation of state < —0.78, 
as implied by Sne la data at 95% c.l. requires either 
a large value for (3, or a very small B. Since the last 
case corresponds to an expansion history of the universe 
practically indistinguishable from GR, any non-standard 
behavior induced by the ST theories in the past should be 
likely accompanied by an equation of state different from 
— 1 today. If DE is a pure cosmological constant, then 
the bound from Cassini implies B < O(10~'^) (making 
A((p) practically indistinguishable from one at least since 
BBN on), or unnaturally large values of 
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FIG. 4: The impact of the Cassini GR test. The regions below 
the curves are excluded at la level. The SNIa bound on the 
DE equation of state is also shown. 



C. CMB power spectrum 



measurements of the peak locations [l9j do not translate 
straightforwardly into a bound on Aifldoc)- For instance, 
it is found that a shift in the peak multipole can be ob- 
tained also by varying the energy densities according to 

m 

AZpeak , oc;^^ 

- ~ -l.Zb—- 

^peak ^ ' 

-0.23 ^ " + 0.09—^ + 0.089-^ 

SO that, in general, a full reanalysis of the CMB including 
the new parameter A(adcc) would be required. However, 
in the present case we find that, once the BBN bound 
on A has been imposed, the resulting values for A(adec) 
are so close to unity as to give shifts in the peak mul- 
tiples smaller than the experimental error. Thus, the 
CMB spectrum does not provide significant bounds to 
the present scenario. 



The impact of a cosmological constant or quintessence 
on the CMB power spectrum has been extensively an- 
alyzed in refs. [2^. In the context of ST theories, the 
problem has been studied in Refs. [TsL l2^. The main 
change with respect to a theory for DE based on GR is 
due to a difference in the expansion rate, which affects 
the angular scale of the anisotropics. The angle under 
which the first peak is seen goes as 

^pcak 7r//poak Vs tdoc Zdcc/d{zdoc) , (27) 

where ^peak is the corresponding multipole, Vg is the 
sound speed, idee and Zdcc are the time and redshift of 
decoupling, and c?(zdcc) the distance to the last scattering 
surface. The latter is given by 

dizdoc)^ I (28) 

Jl/(Zdoc + l) aH 

and is thus dominated by the behavior of H close to 
the upper limit of integration, where aH is smaller. For 
this reason, ST theories passing the GR tests (A ^ 1 
today, that is, H Hq^) imply a small deviation of 
the distance to the last scattering surface with respect to 
GR. 

On the other hand, the decoupling time might be sig- 
nificantly more perturbed. It is given by an expression 
analogous to Eq. H28I) with the upper (lower) limit of inte- 
gration replaced by l/(zdec + 1) (0). As a result, since the 
universe expanded faster than in GR at early times, we 
expect tdcc to be smaller, and the peak to move towards 
higher multiples. We find 

A^pcak 4 A(adcc) - 1 .„„^ 

—, — o — jp — ^ — ' y^^) 

which is consistent with the numerical findings of ref . |l8j | . 

Due to the well known degeneracy of the CMB spec- 
trum with respect to cosmological parameters, present 



V. IMPACT ON WIMP RELIC ABUNDANCE 

Having in mind all the bounds discussed in the previ- 
ous Section, we can now go on to compute the cosmo- 
logical evolution of the scalar field and its impact on the 
DM relic abundance. 

As a first step, we want to estimate if ST can have 
a sizeable effect on the Jordan-frame Hubble parame- 
ter H at the time of WIMP decoupling, without violat- 
ing any of the available cosmological observations. We 
will consider the function A{(p) of Eq. 1)19(1 . imposing 
on the parameters B and /? the phenomenological con- 
straints already discussed. We will then compute the 
ratio H I Hqy\. at the decoupling time of a typical WIMP 
of mass m = 200 GeV. In this way we will be able to get 
an estimate of the effect before going into further detail. 

The tightest bound is that coming from Eq. (|24|) . It has 
an impact on both B in Eq. H19() and on the initial condi- 
tions of ip at temperatures higher than the WIMP freeze- 
out. Indeed, since on the tracker solution the scalar field 
is = 0{1) today, it should have been ^ 1 at nu- 
cleosynthesis, otherwise it would not have reached the 
attractor in time This implies B < 0(0. 1). 

As already discussed, the equation for the dynamics of 
the scalar field tp is obtained by substituting the expres- 
sion of Eq. ((T5|l in the RHS of Eq. ((T^ and choosing a 
coupling function a{(p) as defined in Eq. H20|) . In the sum 
of Eq. 1(1 5|l only the terms corresponding to particles with 
m < Tc have been considered, i.e. particles lighter than 
the critical temperature of the phase transition through 
which they acquire a mass (see Ref. [?!). In particu- 
lar we have taken into account the top quark, the Z", 
the W^, the bottom quark, the tau quark, the charmed 
quark, the pions, the muon, the electron and a WIMP 
particle of mass m — 200 GeV. Numerical integration of 
the equation for (p has been carried out between approx- 
imately 500 GeV and today. We have then computed. 
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FIG. 5: The contours show the expansion rate enhancement 
H/Hgh at T = 10 GeV obtained in the ST model, as a func- 
tion of the initial values of the factor A{ip) and of the ratio 
of the scalar to background energy density pip/pB- We con- 
sidered for the initial conditions a temperature of T = 500 
GeV. The black area represents initial conditions which are 
excluded by nuclesynthesis. The grey contours represent en- 
hancements of 1 10^ 10^ H- 10^ 10^ H- 10*, 10" ^ lif from 
the lightest to the darkest. The dashed lines show the shifts 
of the CMB doppler peaks obtained in the ST model. 



through Ea. H14|l . the modified expansion rate in the ST 
theory at a temperature T = 10 GeV corresponding to 
a typical time of WIMP decoupling and compared it to 
the expansion rate of the standard case at the same tem- 
perature. 

In Fig. iniwe plot the ratio H/Hgr at T = 10 GeV 
as a function of the initial value of A{ip) and initial ratio 
of the scalar to background energy density pip/ pb- We 
have restricted the possible initial conditions to those re- 
gions of parameters values respecting the BBN bound of 
Eq. (|24|l . We see that we have been able to produce an 
enhancement of the expansion rate up to O(IO^) at the 
time of WIMP decoupling. It is then worth studying in 
more detail what happens to the WIMP relic abundance. 

Let us now consider the calculation of the relic abun- 
dance of a DM WIMP with mass m and annihilation 
cross-section {uannv). As already mentioned, laboratory 
clocks and rods measure the "physical" metric g^i, and 
so the standard laws of non-gravitational physics take 
their usual form in units of the interval ds^ . As out- 
lined in Ref.jl2[, the effect of the modified ST gravity 
will enter the computation of particle physics processes 
(like the WIMP relic abundance) through the "physical" 
expansion rate H defined in Eq. . We have therefore 
to implement the standard Boltzmann equation with the 



dY 



1 s 



(31) 



where x ^ m/T, s = (27r2/45) /i*(r) is the entropy 
density and Y = n/s is the WIMP density per comoving 
volume. 

We have considered values of H wich respect all the 
bounds discussed in Section IV. Specifically, we have con- 
sidered the function A{ip) as given in Eq. ifT^ll with pa- 
rameters B = 0.1 and /? = 8. The function A{ip) for this 
choice of parameters is plotted in Fig. El which shows 
that A{ip) is very large at large temperatures, and then, 
at a temperature T^, sharply drops to values close to 1 
before nucleosynthesis sets in. A parametrization of the 
behaviour of A{ip) for T > T^, that will be useful in the 
following discussion, is: 



A{^) 2.19- 10^4 



0.82 



9.65 • 10^ 



GeV\°*^' 



,0.82 



(32) 



where Tq is the current temperature of the Universe. 

We have numerically checked that, in the regime we 
are considering, a good approximation to the physical 
Hubble parameter is given by: 



H - A{ip)HGR 



(33) 



The solution of the Boltzmann equation is therefore for- 
mally the same as in the standard case, with the notice- 
able difference that now the Hubble parameter gets an 
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FIG. 6: A typical behaviour of the function A{ip) defined in 
Eq. 119II . calculated for parameters B = 0.1 and /3 = 8. 
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FIG. 7: Numerical solution of the Boltzmann equation 
Eq. lEU in a ST cosmology for a toy-model of a DM WIMP 
of mass m — 50 GeV and constant annihilation cross-section 
(o"annf) = 1 X 10~^ GeV~^. The temperature evolution of the 
WIMP abundance Y{x) clearly shows that freeze-out is an- 
ticipated, since the expansion rate of the Universe is largely 
enhanced by the presence of the scalar field Lp. At a value 
X = m/T^ a re-annihilation phase occurs and Y{x) drops to 
the present day value. 



additional temperature dependence, given by the func- 
tion A{ip). This can be translated in a change in the 
effective number of degrees of freedom at temperature T: 



7*(a;) — > A'^{x)gi,{x) 



(34) 



An approximated solution of Ea. (|31|l can be cast in a 
form analogous to the standard case: 



1 

Yo 



1 



45G 



m 



1/2 

where G{x) = hi,{x)/gj {x) and Yq and Yf are the 
WIMP abundances per comovin volume today and at 
freeze-out, respectively. The freeze-out temperature is 
obtained by the following implicit equation: 



= In 



0.038Afp5TO 



1727 



1/2 



A{xf)gJ {xf)_ 



(36) 



where g is the internal number of degrees of freedom of 
our WIMP. Clearly, when A{x) — > 1 we recover the stan- 
dard case. The relic abundance is then simply given by: 



Pcrit 



(37) 



where sg is the present entropy density and pcrit denotes 
the critical density. 



A numerical solution of the Boltzmann equation 
Eq. (|31|1 in a ST cosmology with the function A{x) given 
in Fig. El is shown in Fig. [3 for a toy-model of a DM 
WIMP of mass m = 50 GeV and constant annihilation 
cross-section {ca.nn'v) = 1 X 10~^ GeV~^. The temper- 
ature evolution of the WIMP abundance Y{x) clearly 
shows that freeze-out is anticipated, since the expansion 
rate of the Universe is largely enhanced by the presence 
of the scalar field ip. This effect is expected. However, 
we note that a peculiar effect emerges: when the ST 
theory approached GR (a fact which is parametrized by 
A{lp) — ^ 1 at a temperature T^, which in our model is 0.1 
GeV), H rapidly drops below the interaction rate P es- 
tablishing a short period during which the already frozen 
WIMPs are still abundant enough to start a sizeable 
re-annihilation. This post-freeze-out "re-annihilation 
phase" has the effect of reducing the WIMP abundance, 
which nevertheless remains much larger than in the stan- 
dard case. For the specific case shown in Fig.[7|the WIMP 
relic abundance is VlK^ = 0.0027 for GR, while for a ST 
cosmology becomes VLh^ = 0.12, with an increase of a 
factor of 44. 

The phenomenon of re-annihilation can be conve- 
niently discussed in terms of the relation between the 
expansion rate of the Universe H and the WIMP inter- 
action rate P = y s (cTanni') ■ A numerical calculation of 
these two quantities is plotted in Fig. |S1 as a function 
of the temperature. The departure from equilibrium oc- 
curs earlier than in the GR case, because H ^ i?GR- 
When decoupling is completed, the particles evolve with 
an approximately constant Y ^ Yf and P ~ T"^, while 



ol/2 



I.e. 



the Hubble rate evolves as H ~ A{x) p 
slower than in the standard case (we have used here the 
approximate A{x) behavious of Ea. (|32l) '). 

At the transition temperature the Hubble rate 
drops to its standard value i?GR smd becomes smaller 
than the interaction rate: in this case the decoupled 
WIMPs start to annihilate again, for a short period. Af- 
ter this re-annihilation phase, the particles continue to 
evolve with an approximately constant abundace Y < Yf 
and r recovers the behaviour T'^, while Hqr ~ P^^^ ^ 
as usual. 

We notice that a re-annihilation phase does not occur 
in the case of kination, i.e. in the case the energy den- 
sity of the Universe is dominated by the kinetic term of a 
scalar field In this case the evolution of the expan- 
sion rate is -ffkin ^ during kination, and than evolves 
smootly into the standard behaviour i/cR ~ T^- Dur- 
ing kination both H and P have the same T-dependence 
and closely follow each other, until kination ends and the 
standard behaviour is recovered. Re-annihilation is pos- 
sibile if the phase during which the expansion rate has 
the transition toward its standard GR behaviour is faster 
than the post-freeze-out evolution of the interaction rate, 
i.e. faster than T^. 

The change in the freeze-out temperature is shown in 
Fig. |51 where we show the ratio between the freeze-out 
values of Xf = m/Tf in ST cosmology and in GR. The 
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FIG. 8: The Expansion rate of the Universe H and the WIMP 
interaction rate T = Y s {(Tanni') are plotted as a function of 
the temperature. The re-annihilation effect discussed in the 
text is outlined. The small drop in the rates at T = 300 MeV 
is due to the quark-hadron phase transition. 



FIG. 10: Increase in the WIMP relic abundance in ST cos- 
mology with respect to the GR case. The solid curve refers 
to an annihilation cross section constant in temperature, i.e. 
(fannw) = a = 10~^ GeV~^, while the dashed line stands for 
an annihilation cross section which evolves with temperature 
as (o-annu) = b/x = 10"^ GeV'^/x. 
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FIG. 9: The ratio between the freeze-out values of Xf = 
m/Tf in ST cosmology and in GR as a function of the WIMP 
mass. The dashed, solid and dotted lines refer to (aannf ) = 
10"* GeV-^ 10"^ GeV"^ and 10"" GeV-^ respectively. 



between 10% and 40% for WIMPs which can provide 
abundances in the cosmologically acceptable range. 

The amount of increase in the relic abundance which 
is present in ST cosmology is shown in Fig. 1101 The 
solid curve refers to an annihilation cross section con- 
stant in temperature, i.e. (aannv) = a, while the dashed 
line stands for an annihilation cross section which evolves 
with temperature as: (uanni^) = b/x (these two cases 
correspond to the two limiting situations of the usual 
non-relativistic expansion of the thermally averaged an- 



nihilation cross section: (da 



b/x). In the case 



of s-wave annihilation the increase in relic abundance 
ranges from a factor of 10 up to a factor of 400. For a 
pure b/x dependence, the enhancement can be as large 
as 3 orders of magnitude. 

The behaviours shown in Fig. ^| which have been 
obtained by a numerical integration of the Boltzmann 
equation Eq. (|31(l . can be understood by employing the 
approximate analytical solution H35|l . In the case of 

(CTannW) = O, Eq. (jSHJl givCS: 



1 

Yo 



= gmG{xf^) 



(38) 



freeze-out temperature is anticipated about a factor of 2, 
with a dependence also on the annihilation cross section, 
as is clear from Eq. H3t)|) : for very low values of (cTannw) 
the freeze-out temperature may be anticipated up to a 
factor of 5. For these low cross sections the relic abun- 
dance is anyway largely overabundant: we can therefore 
quantify the reduction in in a factor which ranges 



in the standard GR case, and 



1 



— = Qm 



Q m 



G{xf) 



A 1.82 



G{x^ 



(2.ST)1.82 



[x^ 



(39) 
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in our ST model where the A{x) function is given in 
Eq. (EIJ) for T > and A{x) = 1 otherwise {A = 
9.65-103(GeV/m)°-®2). For the sake of simpUcity, in both 
solutions we have dropped the term l/Yj which adds a 
small correction, not relevant for the present approximate 
discussion. In both equations Q = -^/tt/ (45 G). The ratio 
R of the relic abundances is: 



1.82 Ax^ x' 



GR 



1.82 Arc xj 



1.82 



(40) 



x^'^ and we have 



where we have approximated x'^^ 
defined ra = G{x^)/G{x f). By making explicit the mass 
dependencies we obtain: 



R 



Ar TOGcV 

Br + m\,ll 



(41) 



where the mass is expresses in GeV, Ar ~ 1.76 • 
]^04 2,0.82 ^ 2.05 • 10^ Br = 1.76 • lO^T^rcx}-®^ ^ 
2.05 • 10^, and the numerical values have been obtained 
for a;/ ~ 20 and re — 0.5 (since in our case is smaller 
than the quark-hadron phase transition which we have 
set at Tqcd = 300 MeV). The analytic approximation 
of Eq. H41|l helps to explain the behaviour shown by the 
solid curve in Fig. 1101 which has been obtained by nu- 
merical calculations which employ the exact form of the 
function A{ip). From Eq. I|41|) we can in fact derive that, 
for low masses, the ratio R has the behaviour: 



R 



ml ^ Tf 
tgT^ Xf re 



(42) 



which shows that in this mass regime R grows almost 
linearly with the WIMP mass rn, and it is larger for lower 
values of T^. If we accept as low as the BBN scale, 
we can obtain a further increase in the relic abundance 
of a factor 100 on the top of the one showed in Fig. ^] 
for low values of to. When the WIMP mass is very large, 
the ratio R behaves as: 



R 



1.76-104482 



(43) 



with a slight drop with the mass. The position of the 
maximum and the maximal value of R are given by: 



™max 
™GeV 



(2.15-10VgT^)°-56 



and: 



108 



^0.45 



(44) 



(45) 



These expressions show that the maximal effect is also 
obtained for the lowest values of in this case the 



position of i?max is shifted toward lower masses. For 
at the BBN scale, the maximal increase in the WIMP 
relic abundance is of the order of 3000, instead of about 
400 obtained for = 0.1 GeV and shown in Fig. El 
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FIG. 11: Increase in the WIMP relic abundance with an 
annihilation cross section {(Jannf ) = a, for diflerent values of 
a. The dot-dashed, solid and dashed lines correspond to a = 
10"'' GeV-^ 10-'^ GeV-2 and 10"" GeV-^ respectively. 



An interesting property shown by Eq. (|41|l is that R 
does not depend explicitely on the annihilation cross sec- 
tion (cTann^) = which drops out in the ratio. An im- 
plicit dependence on the cross section is present through 
X/, as can be seen in Eq. (|36|l . This dependence however 
is only logarithmic and does not spoil the general be- 
haviour of R shown in Fig. ^| This is shown in Fig. ITTI 
the largest difference occurs for very low annihilation 
cross sections, for which the deviation of xj is larger. 
However, for cross sections of interest, i.e. cross section 
which provide relic abundances below the cosmologically 
acceptable upper bound, the values of R are stable to a 
relatively good extent. 

A similar analysis holds in the case of (cTannf) = b/x. 
However, in this case the dependence of R with Xf is 
somehow stronger (as obtained from the integration in 
Eq. ||23)), and the effect of changing ((Tann'f^) is slightely 
larger. This effect can be seen in Fig. El where R is 
shown for different values of the parameter b. Notice that 
larger cross sections, which in the standard case provide 
lower values for the relic abundance, are the ones which 
get more enhanced in ST cosmology. 

Finally, as an example we show in Fig. El the relic 
abundance as a function of the WIMP mass in the case 

of (cTannW) = a — 1 ■ 10^^ GcV and (iTannV) = b/x = 

1 ■ 10^** GeV/x. We see that, in this case, the relic abun- 
dance can be at the level required to explain the CDM 
content of the Universe (fJcDM^^ = 0.095-^-0.13 2M) ^r 
a ST theory, while it is underabundant in the standard 
case. The models shown in Fig. El represent a case in 
which we can explain at the same time both the DM and 
DE contents of the Universe, and the interplay of the two 
component is crucial in determinig the right abundances 
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FIG. 12: Increase in the WIMP relic abundance with an 
annihilation cross section (crann^') ~ b/x, for different values 
of b. The dot-dashed, solid and dashed lines correspond to 
b = 10-* GeV-2 10"^ GeV-2 and 10"^° GeV-^ respectively. 
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FIG. 13: The relic abundance in a ST theory as a function 
of the WIMP mass in the case of ((Tannw) = a = 1 • 10~^ 
GeV (solid line) and {a ) =b/x^l- IQ-'^ GeV/x (dashed 
line). The upper (lower) dotted lines corresponds to the GR 
case for (aannv) = a — 1 ■ 10~^ GeV and (cTannw) = b/x — 
1 • 10~* GeY/x, respectively. 



of both DM and DE. 

An analysis of specific particle candidates of DM, in 
particular in supersymmetric models, will be examined 
elsewhere 1221. 



VI. CONCLUSIONS 

The idea of exploiting primordial (ultralight) scalars 
in order to shed some light on a dynamical interpre- 
tation of DE is by now a widespread research topic in 
the literature. In this paper we follow the promising 
proposal of considering the quintessence scalar as em- 
bedded in a scalar-tensor theory of gravity. This ap- 
proach is at variance with the usual interpretation of 
quintessence as a new light scalar whose interactions with 
matter are subject to the tight phenomenological con- 
straints on the equivalence principle violation and time- 
variation of the fundamental coupling constants. Iden- 
tifying the quintessence field with the scalar component 
of a ST theory, instead, does not pose any threat on the 
equivalence principle, since by construction matter has a 
purely metric coupling with gravity. 

We focus on quintessence ST models which possess a 
double "attraction mechanism" , one to GR and the other 
ensuring puE to follow a tracking solution. These two 
simultaneous mechanisms act as a "protection" for the 
theory to prevent its fall into immediate troubles (for 
instance, large departures from GR predictions). Nev- 
ertheless, we still obtain important phnomenological sig- 
natures which might disentangle this theory from GR or 
other alternative proposals. The tests of our ST scenario 
divide into two classes: deviations from GR and depar- 
tures from standard cosmology, in particular concerning 
the expansion rate of the Universe. 

The latter effect has a big impact on the most distant 
epoch of the Universe for which we have "direct" informa- 
tion, i.e. nucleosynthesis. However, we pointed out that 
we can further extend the implications of a non-standard 
H in the early Universe to times prior to nucleosynthe- 
sis. Sticking to the standard WIMP picture of DM, one 
of the most relevant events before BBN is the WIMP 
decoupling which is expected to have occurred at a tem- 
perature of a few GeVs. Our work shows that, despite 
the severe "filters" on ST quintessence models which are 
provided by BBN, solar system tests of gravity and, to a 
lesser degree, by GMB, it is still possible to find remark- 
able enhancements on the expansion rate of the Universe 
at WIMP freeze-out, yielding to relic WIMP abundances 
which can vary up to a few orders of magnitude with 
respect to the standard case. 

In this paper we pointed out some general features of 
the new "WIMP story" around its decoupling tempera- 
ture in the presence of ST quintessence. In particular, we 
noticed that some unexpected effect can take place, such 
as a short phase of WIMP "re-annihilation" when ST 
approaches GR. Needless to say, such potentially (very) 
large deviations entail new prospects on the WIMP char- 
acterization both for the choice of the CDM candidates 
and for their direct and indirect detection probes. A thor- 
ough reconsideration of the "traditional" WIMP identi- 
fied with the lightest neutralino in SUSY extensions of 
the SM as well as the identification of other potentially 
viable CDM candidates in the ST context is presently 
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